Nuclear a-particle condensates: Definitions, occurrence conditions, and consequences 
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There has been a recent flurry of interest in the possibility of condensates of a-particles in nuclei. 
In this letter we discuss occurrence conditions for such states. Using the quantality condition of 
Mottelson we show that condensates are only marginally expected in a-particle states. We proceed to 
demonstrate that few-body nuclear condensates are ill-defined, and emphasize the conflict between 
Q-localization and a-condensate formation. We also explore the connection between Ikeda diagrams, 
linear chains, and Tonks-Girardeau gases. Our findings show that no new information is contained 
in the approximations of nuclear states as a-cluster condensates. Furthermore, condensates of more 
than three a-particles are very unlikely to exist due to couplings to other degrees of freedom. 

PACS numbers: 03.75.Hh, 21.45.+V, 21.60.Gx 



Introduction. The idea of a-particles as essential con- 
stituents in the structure of nuclei arises from the small 
radius, the relatively large binding and the spin satura- 
tion of both neutrons and protons. Attempts were made 
in the early days of nuclear physics to construct nuclear 
structure from a-particles and valence nucleons [jj, [2], l3j • 
In general these attempts were largely unsuccessful be- 
cause the nucleon distances within and between different 
a-particles are comparable. Thus there are no compelling 
reasons for clusterization of nucleons into a-particles. 

On the other hand, a-cluster models were able to ex- 
plain many properties of specific light nuclei This 
was high-lighted by the prediction of a 3a-structure near 
threshold [5| which is crucial for the nuclear synthesis of 
the heavy elements in stars. This Hoyle state was soon 
found experimentally [|| and its properties established in 
a microscopic cluster model [3] where it was first charac- 
terized as a "gas-like" structure. This structure is con- 
firmed in details in numerous theoretical works 0, [Io| . 
A radius about 30% larger than that of the ground state 
of 12 C is also reproduced by use of an approximate wave 
function consisting of an antisymmetrized product of 
three gaussians each containing an a-particle [ill ]. Also 
the ground state of 8 Be was in the same approximation 
described as a gas- like structure of two a-particles [T^ |. 
and the established structure [H, [T3| was again essen- 
tially recovered. The fundamental continuum resonance 
properties were suppressed by use of box-like boundary 
conditions effectively supplying a confining external field. 

These approximations were presented as novel dis- 
coveries of condensates consisting of two and three a- 
particles. The last five years have witnessed surprisingly 
large efforts invested in both investigations of the accu- 
racy of the approximation in [111 ] and extensions to simi- 
lar simple models for other nuclei (l5l[T3|. The aim seems 
to be a search for a-cluster condensates in nuclei. The 
inspiration is from atomic physics where Bose-Einstein 
condensates (BEC) of cold atoms are routinely made and 
manipulated by external fields [l7|. Related theoretical 
investigations are also abundant, see e.g. the review [l8| . 

The concept of BEC is well-defined for macroscopic 
systems of cold atoms and molecules. Extensions to self- 



bound quantum systems with a small number of parti- 
cles are not straightforward. The purpose of this let- 
ter is to discuss the concept of few-body nuclear con- 
densates, give definitions, compare to cold atomic gases, 
show the conflict between localized a-cluster models and 
a-condensation, derive occurrence conditions and inves- 
tigate consequences. 

Concepts. To find a-cluster condensates in nuclei two 
conditions must be met, i.e. (i) the nucleons must be con- 
fined in a-clusters, and (ii) these a-particles must form 
a condensate. Although a-cluster models in general are 
unsuccessful we shall assume that (i) holds. Then to en- 
sure nucleon antisymmetry the a — a distance must on 
average be larger than the diameter of the a-particle. 

The classical definition of an ideal Bose-Einstein con- 
densate is a collection of identical particles in the same 
quantum state. This implies that the independent parti- 
cle model gives an accurate description which obviously is 
true for non-interacting particles in an external field. For 
interacting particles the mean-field description is valid 
when Mottelsons quantality condition is met pj| . i.e. 
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where m is the mass of the particles and c m i n is the 
distance between two particles when the total two-body 
potential has its minimum value V m i n . When Ajvf < is 
small the attractive potential dominates over the kinetic 
energy and the particles are confined to the attractive 
pockets, i.e. localization or solid structure. When Ajv/ot 
is large the particles can not be confined by the attraction 
and the mean-field model is appropriate. 

The condition for Bose-Einstein condensates is that the 
deBroglie wavelength A^s of the motion must be larger 
than the distance c m i n to the nearest neighbor, i.e. 



1 < A bec EE 



A 



dB 



2irh 



= 7TV2A 



Mot 



(2) 



Thus A^ ec = 2ir 2 KMot implying that the quantality in- 
equality in eq.([T]) separating solid and mean-field struc- 
tures is equivalent to the condition for breakdown of the 
classical gas regime in statistical mechanics. 
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We can evaluate these conditions for the a — a poten- 
tial V without any bound states parametrized in [201 ] as 
an attractive and a repulsive gaussian of different ranges 
(r a ,r r ) and strengths (V a ,V r ), i.e. 



V(r) = V r exp(-r 2 /r 2 ) - Kexp(-r 2 /r 2 ) , 



(3) 



where the minimum V m in ~ 5 — 8 MeV (including the 
Coulomb energy of w 2 MeV) for c m i n « 2.5 — 3.0 fm 
and h-Mot ~ 0.1 — 0.2. Thus a-particles would be in the 
mean-field range but with a strong tendency to localize. 

Symmetry requirements. Wave functions describing 
self-bound few-body structures must be invariant under 
translations and rotations. The connection to conditions 
for condensate formation is most easily illustrated by use 
of an iV-body wave function '5 expressed as products of 
identical gaussian single-particle wave functions, i.e. 

N 

*({r,}) = (bV^r 3N/2 exp(- ]T r 2 /(26 2 )) , (4) 

i=i 

where is the i'th coordinate. This as well as all 
other mean-field wave functions violate translation in- 
variance, or equivalently momentum conservation, which 
is restored by integrating ^(jr, — J?'}) exp(iP ■ Ft') over 
all R' . The solution ^>i n t of lowest energy has P = 
which for eq.© results in 

* int ({n}) = (bV^r 3{N - 1)/2 exp( V/(26 2 )) , (5) 



p 2 -E^ = ^E4=E^ 2 -^ (6) 
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qi = ri R , rij=ri-Tj , R = — n , (7) 



where the coordinates now are measured from the com- 
mon center-of-mass (cm) R. This wave function is in- 
variant under rotations around the cm. 

In contrast to the mean-field solution the particles can 
be correlated and the wave function ^i oc in a body-fixed 
coordinate system localized in distributions around pre- 
ferred points .Rfc, i.e. 

N 

*toc({ri}) « J2 ex P(- E< r * - R P{1 )) 2 /(2B 2 )) , (8) 
p 
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where the normalization is omitted and full symmetry is 
achieved by the summation over all permutations p of the 
set of numbers {1,2,...., N}. The translational invariance 
is restored precisely as for eq.©, i.e. the wave function is 
obtained from eq. © by the substitution — > qi in eq. © 
and a corresponding change of normalization constant. 
The rotational invariance is broken for ^i oc in eq. © but 
recovered for states of zero angular momentum by linear 
combinations of all spatial rotations of ^i oc . 



Condensate assessment. To decide if a given wave 
function describes a condensate we apply different avail- 
able definitions. We illustrate again with gaussian wave 
functions which overestimates the degree of factorization 
of the JV-body wave function. One necessary criterion for 
a condensate is that the one-body density matrix must 
have an eigenvalue A comparable in size to N [25j. For 
a mean-field product solution the condensate fraction is 
Cf = X/N = 1. However, the one-body density matrix 
is ill-defined for self-bound systems of a finite number of 
particles. This is due to the cm- motion which decouples 
completely for correct translationally invariant solutions. 
An appropriate cm wave-function could be chosen to al- 
low the usual definition of the one-body density matrix. 
The choice could be such that the condensate fraction 
Cf is optimized which would be equivalent to adding an 
external field as in atomic physics. For eq.([5]) this recov- 
ers the product wave function of all coordinates in eq. ^ 
where the cm-motion is completely ignored. 

Instead of using all particle coordinates relative coor- 
dinates could be used and an internal one-body de nsity 
matrix, n(q,q'), defined 0, M, HI HI • Following [if, 
i.e. inserting q N = - Y^iLi 1 1i m c q-©> we get 

n(q,q')(X J d :i q 2 d 3 q 3 ...d 3 q N _ 1 \^ int \ 2 
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where q and q' refer to q±. The condensate fraction 
obtained through the largest eigenvalue is then [24| c/ = 
8/(1 + y/2 - 2/Nf which decreases with N from 1 for 
N = 2 towards about 0.57 for large N. However, the 
choice of internal coordinates is arbitrary [24[ and we 
could as well choose q\ supplemented by a set of N — 1 
independent Jacobi coordinates. Then the density matrix 
corresponding to eq.© would factorize and give c/ = 1. 

For a-clusters these options can by appropriate choices 
lead to large condensate fractions for rather accurate 
cluster wave functions. This is because approximate fac- 
torization easily arises at smaller distances where the po- 
tential minimum resembles a harmonic oscillator and the 
related s-wave solutions resemble gaussians. 

Instead of using the eigenvalues of the density matrix 
a condensate criterion could be that the one-body (in- 
ternal) density matrix should factorize at large distances 
[25| . This criterion is extremely difficult to fulfill because 
the correct nuclear wave functions never factorize at large 
distances as shown in (26J. Thus, at best only properties 
at intermediate distances could possess condensate prop- 
erties with this criterion. 

Yet another condensate criterion is that all particles 
occupy the same quantum state (24[. This implies that 
removal of one particle should leave the single-particle 
wave functions completely unchanged. However, for a 
finite number of particles the remaining interacting par- 
ticles would reorganize into a different structure. This 
criterion would be extremely difficult to test directly. 
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Localization. The a-cluster models and the quan- 
tality parameter in eqs.([T]) and @ suggest that local- 
ization, or crystal features of the wave function, may be 
important. The resulting condensate fraction depends 
strongly on the degree of localization as we can see ex- 
plicitly by computing the one-body density matrix for 
eq.®. We assume very narrow non-overlapping gaus- 
sians and an appropriate cm-motion, and obtain 



i{r,r') = {B^) 
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(10) 
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x exp(-((r - R k f + (r' - R k ) 2 )/(2B 2 )) , 

k=l 

which has N equally large eigenvalues while all others are 
zero. This is a condensate fraction of 1/7V corresponding 
to one single-particle state for each of the N particles. 
However, after restoration of rotational symmetry only 
eigenvalues zero remain. If the widths, B, of the gaus- 
sians increase and they begin to overlap with each other 
one eigenvalue separates out and becomes finite. Increas- 
ing the width leads to increasing overlap with a product 
wave function like eq.Q. We can quantify by comput- 
ing the overlap between the factorized and localized wave 
functions in eqs.(|3]) and ([5]), i.e. 
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which only is close to unity when b ~ B and cither 
Ri/B -C 1 or Ri/b < 1. Eq.(|TTJ| is also obtained by 
replacing \&; oc with the rotationally invariant wave func- 
tion. Thus a substantial condensate fraction requires 
that the overlap with eq.® is large. However, the spa- 
tial extension must be large to ensure definition (i) of 
non-overlapping a-particles. 

Condensates from cluster models. The well-known 
structure of the Hoyle state in 12 C has about 90% overlap 
with a-particles in relative s-waves around the center- 
of-mass [HI, Hi], [22], J27J • This corresponds to an eigen- 
value of about 0.7 |2ll . |22| in agreement with our up- 
per bound of cr = 0.80 derived from eq.©. At the 
same time a-cluster models show a-particle density dis- 
tributions localized around specific points in space [27| . 
Reconciling these results, where apparently both local- 
ization and large condensate fraction are present in the 
same wave function, is only possible with large widths 
of the localized wave in eq.(8). This effectively recovers 
the independent particle wave function in eq.(4) where 
the a-particles are sufficiently separated to remove the 
need for nucleon antisymmetrization. These arguments 
show that the a-condensate states proposed should be 
regarded as merely an approximation to existing nuclear 
a-cluster states. 

A crucial question is whether a condensate structure 
can be experimentally distinguished from other struc- 
tures. To address this question computed and measured 
electron scattering on 12 C was compared in |27| . The 
conclusions are that a-cluster models of the Brink-type 



[ij and the a-condensate states of 11 1 predict virtually 
indistinguishable cross sections and charge distributions. 
In addition, these models and results from more elabo- 
rate microscopic calculations [27j give precisely the same 
charge density at large distances. Thus the classical clus- 
ter parametrization supplemented by nucleon antisym- 
metry [ll| is apparently accurate to about 90% for the 
Hoyle state. However, it is important to realize that ful- 
filling an ambiguous definition has very little to do with 
true condensates which only can be diagnosed through 
properties of the wave functions and not by density dis- 
tributions. In particular, observable coherence properties 
of the many-body states are necessary to separate clus- 
ter states from condensates. Both this and the localiza- 
tion discussion above strongly indicate that no new con- 
sequences arise from approximating cluster states with 
"condensates" . 

Condensate identity. The approximation as a con- 
densate wavefunction of a quantum state rapidly gets 
invalid with increasing nuclear mass. This is seen from 
a sequence of four arguments. First the approximation 
as a condensate wavefunction is related to a restricted 
part of the full Hilbert space. Variational computations 
of condensates assume a class of wave functions with pa- 
rameters determined by minimizing the energy. When 
the Hilbert space is extended to include other degrees of 
freedom the solution must remain essentially unchanged. 
An analogy is found in the s-wave neutron strength func- 
tion which is broad and distributed over a large number of 
many-body states. This is reflected in the lack of neutron 
halo states at excitations around the neutron binding en- 
ergy B n [H, [29[. To maintain the condensate character, 
the residual coupling V Cj „ of an a-condensate state |c) to 
the true many-body continuum nuclear states \n) must 
all be smaller than their energy difference [29j |. 

The second argument in the sequence is that this ap- 
proximation gets increasingly worse with increasing ex- 
citation energy because the density of states increases. 
A "clean" condensate wavefunction must be more and 
more "smeared out" over the true many-body states and 
at some excitation energy the condensate wavefunction 
no longer describes a state of the nucleus. Third, alpha- 
condensates are postulated at the threshold for disinte- 
gration into alpha-particles. Fourth, this is at an exci- 
tation energy E* of about 7 MeV for 12 C and increas- 
ing by about 7 MeV for each additional a-particle, i.e. 
E* w 7(^4/4 — 2) MeV. 

We first estimate an average V av of V Ct „ by using a 
nucleon-nucleon potential of range b and strength Vq 
between the two states with similar radii R = bA 1 / 3 . 
Both residual kinetic and potential energy contributions 
are then proportional to the number of nucleons A, i.e. 
\V av \ ~ ASq. The radii in \c) and |n) must be comparable 
if the attraction of short range has to keep the conden- 
sate spatially confined in competition with the repulsive 
Coulomb interaction. The condition for maintaining the 
condensate character is then \V av \ < D, where D is the 
average level spacing. 
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We also estimate V av by replacing \c) by the state 
|a, (A — 4)), consisting of an a-particle and the ground 
state, | (A — 4)), of the (A — 4)-system. These wave- 
functions are similar, because the condensate consists of 
alpha-particles, but they are clearly not identical, since 
the ground state of A-4 cannot be accurately described 
by an alpha-cluster model. However, they are similar in 
the approximations of harmonic oscillator potentials or 
gaussian wavefunctions. The differences in the coupling 
matrix elements from using |c) and \a, (A — 4)) can then 
on average be expected to deviate much less than an or- 
der of magnitude. With many states \n) in the average 
this leads as in [2{J [3(| to the estimate of V av ~ W a 
where W a is the strength of the i mag inary part of the 
a-nucleus optical potential. As in [29| we conclude that 
the condensate resembles one of the many-body states 
when W a < D. 

From the imaginary a-nucleus potential we can esti- 
mate the spreading width T a of an a-particle state on 
the true many-nucleon states. For nucleons the spread- 
ing width of a single-particle state of energy e, r sp , is 
in Fermi-liquid transport theory [3l|, [H, [H[ found to be 
T sp = ((e - /i) 2 + 7r 2 T 2 )/r where fi is the Fermi level 
and T the temperature. The constant Tq is related to 
the imaginary potential and the estimate w 33 MeV re- 
sults in r sp « 1.5 MeV for T = and an energy equal to 
the nucleon separation energy [32l . |33| . Analogously we 
estimate T a for a-particles moving in the medium of nu- 
cleons. The a-nucleus separation energy is about 7 MeV, 
which with Tq w 33 MeV again results in r a = 1.5 MeV 
for T — 0. The finite temperature is obtained from the 
average excitation energy, E* — aT^, at the threshold 
for fragmentation into free a-particles. With the level 
density parameter a — ^4/10MeV we get T c w 4 MeV 
y/l — 8/A where A is the nucleon number. In total we get 
the estimate r a «6 MeV almost independent of nucleon 
number. This correponds to a strength W a ~ 3 MeV for 
the appropriate increasing excitation energy. 

We estimate D in the Fermi gas model adjusted phe- 
nomenologically to excitation energies E* rs B n , i.e. 
Dm D cxp(-2y/a(E* - 2A)), where A 12 MeV/VJ 
is the pairing gap A. The level spacing, Do ss 20 MeV /A, 
for E* = 2A is essentially equal to the single-particle 
level spacing which is the correct limit. The extremely 
simple expression for D can only be an average over many 
nuclei at energies where many excited states are present. 

The conditions \V av \ < D and W a < D are then 

S AoiW a < D cxp(-2^/a(E* - 2A)) . (12) 

A very low limit of both So and W a is 1 MeV [H]. With a 
very small value of E* = 2A we get the conservative esti- 
mate of preservation of condensate identity A < \/20 < 5 
or A < 20, respectively. The spreading width estimate, 
even reduced by a factor of 2, is also larger than the level 
distance for A < 14. These estimates are valid when a 
sufficient number of excited states contributes in the av- 
erage around the threshold energy. This is fulfilled for all 
nuclei heavier than 12 C, including 16 0. These conditions 



for survival of the condensate structures are almost al- 
ways violated, and the violation increases exponentially 
with excitation energy. 

Tonks- Giradeau structures. The linear chain struc- 
tures of a-particles at the break-up threshold, Ikeda 
diagrams [351 ] . are conceptually similar to the one- 
dimensional atomic condensates called Tonks-Giradeau 
structures [3(1 l37l ]. The latter have been realized with 
bosonic Rubidium atoms in optical traps that have strong 
repulsive interactions in the ID geometry. This near im- 
penetrability makes the system behave like a ID Fermi 
gas in many aspects. This is analoguous to impenetrable 
a-particles in ID cluster structures. The corresponding- 
states have been searched for and for many years the 
Hoyle state in 12 C was the favorite candidate. This state 
is now claimed as a condensate with a completely differ- 
ent structure. No observable has been found to distin- 
guish between these structures which in any case both 
are approximated as three-body cluster states. 

To assess whether such linear structures could exist in 
nuclei we turn to the two conditions in eqs.([T]) and 
The a — a Coulomb energy is unimportant compared 
to V in eq.Q. It does not change the condition but it 
also cannot provide the confining external field allowing a 
mean-field condensate-like solution. Hence a linear chain 
structure is only possible with localized a-particles. The 
linear chain structure may also be destroyed by couplings 
to other degrees of freedom. As for three-dimensional 
condensates we estimate the survival probability to be 
very small for any excitation energy above 2A. 

Conclusions. The existence of Bose-Einstein conden- 
sates of a-particles assumes first that nucleons clusterize 
into a-particles, and second that a condensate is formed. 
The quantality condition of Mottelson indicates that a- 
particles marginally prefer independent particle motion 
over correlation. We show that definitions of conden- 
sates of very few particles are ambiguous and lead to 
disparate condensate fractions. The origin is conflicts 
between mean- field solutions, correlations and transla- 
tional or rotational symmetries, and between definitions 
related to short and long-range behavior. The differences 
between nuclear and atomic condensates are few versus 
macroscopic number of particles, dilute versus high den- 
sity, self-bound system versus external confining field, 
and ambiguous versus rigorous definitions. 

In conclusion, we have found that the concept of a 
nuclear condensate is of little use. The recent theoreti- 
cal claims of nuclear a-condensates refer to well-known 
cluster states and can be regarded as merely an approx- 
imation to such states. No observable differences can 
be constructed to distinguish these alleged novel struc- 
tures from ordinary cluster states. Occurrence of one 
and three-dimensional nuclear condensates of more than 
3 particles at higher excitation energies are very unlikely. 
They would either be completely unstable and vanished 
into the continuum, or the a-condensate structure would 
cease to exist due to spreading over many-nucleon states. 
In any case if traces remain they are nothing else than 
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parts of ordinary en-cluster states. 
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